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ABSTRACT 

In this paper I analyze the process of formation of thin current structures in the mag- 
netosphere of a conducting accretion disk in response to the field-hne twisting brought 
about by the rotation of the disk relative to the central star. I consider an axisymmetric 
force-free magnetically-linked star-disk configuration and investigate the expansion of 
the poloidal field lines and partial field-line opening caused by the differential rotation 
between the star and a nonuniformly-rotating disk. I present a simple analytical model 
that describes the asymptotic behavior of the field in the strong-expansion limit. I 
demonstrate the existence of a finite (of order one radian) critical twist angle, beyond 
which the poloidal field starts inflating very rapidly. If the relative star-disk twist is en- 
hanced locally, in some finite part of the disk (which may be the case for a Keplerian disk 
that extends inward significantly closer to the central star than the corotation radius), 
then, as the twist is increased by a finite amount, the field approaches a partially-open 
configuration, with some field lines going out to infinity. Simultaneous with this par- 
tial field opening, a very thin, radially extended current layer forms, thus laying out a 
way towards reconnection in the disk magnetosphere. Reconnection, in turn, leads to a 
very interesting scenario for a quasi-periodic behavior of magnetically-linked star-disk 
systems with successive cycles of field inflation, opening, and reconnection. 

Subject headings: accretion, accretion disks — magnetic fields — MHD — stars: mag- 
netic fields 



1. Introduction 

Magnetic processes taking place in the magnetosphere above a thin accretion disk play an important 
role in establishing the structure of disk outflows (winds and jets), in regulating the accretion flow, 
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Fig. 1. — Axisymmetric magnetically-linked star-disk system. 

and in angular momentum transfer. Of particular interest is the situation where there is a direct 
magnetic connection between the star and the disk (the so-called magnetically-linked star-disk 
system, see Fig. 1). This connection may lead to a direct angular momentum exchange between 
the disk and the star and is thus relevant for neutron star spin-up/spin-down events (Ghosh Sz 
Lamb 1978, 1979; Wang 1987; Lovelace et al. 1995). It is also important because it provides a 
mechanism for direct channeling of accretion flow onto the polar regions of the star (Bertout, Basri, 
& Bouvier 1988; Lamb 1989; Patterson 1994; Konigl 1991). 

The evolution and even the sole existence of such a configuration depends critically on the con- 
ducting properties of the disk, the star, and the overlying magnetosphere. This complex physical 
problem can be studied on various time scales; one should start, however, with the shortest rel- 
evant time scale, namely, the rotation period. On this time scale both the central star and the 
disk can usually be approximated by ideal conductors (with the exception of the case when the 
central object is a black hole), and so can the low-density magnetosphere (or corona) that lies 
above the disk. In this case there is no steady state because the differential rotation between the 
disk and the star leads to continuous twisting of the field lines. The magnetic field in the magne- 
tosphere responds to this differential rotation by rapid expansion driven by the increased toroidal 
field pressure. The field lines become elongated along a direction making a roughly 60° angle with 
the rotation axis. This expansion process has been studied in some detail in the framework of the 
force-free model.^ In particular, semi-analytic self-similar models (van Ballegooijen 1994, hereafter 



^The validity of the force-free approximation in these systems is justified by the fact that, due to very low density 
in the magnetosphere, the Alfven velocity there greatly exceeds both the sound speed and the rotation speed. 
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Fig. 2. — Field- line expansion caused by the relative rotation of the disk with respect to the star in 
the case of the n = 0.5 self-similar model, where ^'disk '"^ f~"' (This figure is taken from Uzdensky 
et al. 2002a). 



VB94; Lynden-Bell & Boily 1994, hereafter LBB94; Uzdensky et al. 2002a) have shown that the 
field-line expansion leads to an effective opening of the field lines after a finite (a few radians) twist 
angle (sec Fig. 2). This finite-time singularity has also been observed in numerically constructed 
non-self-similar sequences of force-free equilibria (Uzdensky et al. 2002a) and also in non-forcc-free 
full-MHD numerical simulations by Goodson et al. (1997, 1999). This process is essentially very 
similar to the process of finite-time field-line opening of a twisted coronal field studied extensively in 
solar physics (see, e.g., Barnes &; Sturrock 1972; Low 1986; Roumeliotis et al. 1994; Mikic &; Linker 
1994; Wolfson 1995; Aly 1995). Thus, at present there is growing evidence that a finite-time open- 
ing of field lines in response to differential rotation is a generic feature of force-free magnetospheres 
of magnetically-linked star-disk systems. 

It is, however, still not clear what happens after this opening. Basically, one can envision two 
possible outcomes: the field lines may either reconnect (across the separatrix between the two 
domains with the opposite direction of the radial field) and close back; or they may stay open 
indefinitely. Li the first scenario (Aly & Kuijpers 1990, VB94, Goodson et al. 1997, 1999) some 
resistivity is present in the magnetosphere, and as a result, a small amount of flux gets immediately 
reconnected, leading to the formation of an X-point. As reconnection progresses, a toroidal plasmoid 
forms in the magnetosphere. This plasmoid contains most of the toroidal flux; it becomes completely 
detached from both the disk or the star and just floats away, presumably feeding a jet (Goodson 
et al. 1999). At the same time, the magnetic link between the star and the disk is reestablished, 
with the reconnected magnetic field lines being now twisted less than they were just before the 
reconnection event. The poloidal magnetic field tension then quickly contracts them back towards 
the star, and the system returns to its initial state. Thus, in this scenario the situation is manifestly 
time-dependent; the time evolution consists of periodic cycles of successive field line infiation due 
to the differential shearing, "effective" field-line opening, reconnection, and the relaxation of the 
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reconnected magnetic field to the initial state (of lesser magnetic stress), which completes the cycle. 

In the second scenario, investigated by Lovelace et al. (1995), there is no reconnection and a true 
steady state is established, whereby the field is at least partially open and the magnetic link between 
the star and a significant part of the disk is permanently broken. 

The fact that the above two scenarios differ so dramatically raises the level of urgency of identifying 
the conditions under which each of them can be realized. The main difference between these two 
pictures is that magnetic reconnection is allowed to take place in one but not in the other. Thus, 
as it often happens in astrophysics and space physics, the physics of magnetic reconnection plays a 
key role. This issue is by no means trivial and it is questionable whether it could be resolved even 
by a direct numerical simulation. Indeed, numerical simulations typically suffer from unrealistically 
large numerical resistivities that make it very easy for oppositely directed field lines to reconnect 
(as seems to be the case with the simulations by Goodson et al. 1999). 

It is usually the case that in order to get efficient reconnection, a thin current sheet is needed. In the 
two-dimensional axisymmetric situation discussed here, a natural place where such a current sheet 
can arise as the system approaches the open-field state is the current concentration region between 
the two domains with oppositely directed radial magnetic field. While the field lines are still closed, 
the current is concentrated along the apex Hue 9 = ^ap(*)-^ As the field opens up, this apex line 
turns itself into the separatrix between two regions of oppositely directed, open (i.e., extending 
to infinity) field lines that comprise the open-field configuration. If one now considers such a 
configuration (where the field lines have already been opened), then one discovers that the formation 
of a current sheet is essentially unavoidable, at least as long as the force-free approximation is valid.^ 
Indeed, an open magnetic field is potential, it has no toroidal field. This is because all toroidal 
flux created by the initial twisting has escaped to infinity as the result of magnetic field expansion 
and opening. The field becomes predominantly radial everywhere, with Br reversing across the 
separatrix 9 = ^ap- In the absence of the toroidal field at 9 = 6'ap, the pressure balance across 
the separatrix cannot be maintained, and this leads to the collapse of the magnetic configuration 
to one with an infinitesimally thin current sheet; the two oppositely directed magnetic fields move 
toward each other, forming a thin current layer. Finally, the collapse is stopped when resistive (or 
other non-ideal) effects become important in the layer. Because the current density is tremendously 
increased, magnetic reconnection can occur. 



^Here, ^ap(^) marks the angular position of the apex of an inflated fleld hne (an apex is deflned as the farthest 
from the central star point on a field line). Usually, 9a,p is close to 60° and depends only very weakly on the field 
hne 

^As was pointed out by the referee of this paper, if a heated atmosphere is present above the disk, then, once the 
magnetic field opens up, the drop of magnetic intensity with distance may be fast enough for the disk wind far out 
to dominate and keep the field from reconnecting. The centrifugal wind could also hold the field open provided there 
is a lower, sufficiently heated atmosphere to feed that wind. This may be the implied reasoning behind the apparent 
neglect of the possibility of reconnection by Lovelace's et al. (1995). 
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It is important to realize, however, that, strictly speaking, as long as one has a force-free magnetic 
field with the field lines closed, one can never have a true, infinitesimally thin current sheet. Indeed, 
while the system is going through a sequence of equilibria of closed field lines, i.e., for t < tc (where 
tc marks the moment time of field-line opening), there is always finite toroidal field, 5^, present. 
This is because it is the toroidal flux that drives the expansion along 6 = O^p, and in a force- free 
equilibrium its outward pressure balances against the poloidal field's large curvature. Thus, in a 
closed force-free configuration the toroidal flux is kept in place by the poloidal field tension and 
prevents poloidal field lines from contracting back to a less-stressed state. At the same time, this 
toroidal field provides the pressure in the 6 direction, which prevents the poloidal field collapse into 
a current sheet configuration. 

One could ask, however, whether a current sheet can form asymptotically, that is, whether the 
sequence of force-free equilibria, which governs the system's evolution, can asymptotically lead to 
stronger and stronger thinning of the current concentration region as one approaches the critical 
moment, so that the characteristic angular width of this region, A6, goes to zero as t ^ tc- Note 
that this does not always happen. For example, in the self-similar force-free model for a uniformly 
rotating disk (VB94; Uzdensky ct al. 2002a), where all magnetic quantities are power laws of r 
with fixed power-law indices, does not go to zero but approaches a finite value that depends 
solely on the flux distribution ^'disk(^) on the disk surface. In particular, if ^'disk('') ^ r^", then 
A9 approaches a finite value of order a fraction of one radian if n = 0(1) (and is proportional to n 
in the limit n — ^ 0), as the critical moment tc is approached. 

Asymptotic current sheet formation was in fact observed in the self-similar model by LBB94 (in 
cylindrical geometry) and also in an essentially very similar work by Wolfson (1995) (in spherical 
geometry in the solar corona context). This is explained by the fact that the power law index n was 
forced to change during the evolution. In particular, in the cylindrical case considered by LBB94, 
all magnetic flux was required to go through the boundary at some small radius (compared with 
the radius under consideration). Therefore, for each value of the power exponent n, the resulting 
equilibrium corresponded to a member of the VB94 family of solutions in the limit of infinite 
expansion, t — > tc{n) (this limit was analyzed in detail by Uzdensky et al. 2002a). Since in the 
VB94 analysis tc{n) grows as n ^ 0, the gradual increase of the twist angle with time meant that, 
when considering the LBB94 solutions, one were bound to obtain a sequence of solutions with 
ever diminishing values of n [determined from the condition tc{n) = t]. That is, in the LBB94 
sequence of solutions, n had to decrease and approached zero at some point, which lead to the 
reported thinning and formation of a current sheet. Related to this change of n was the fact that 
the field-line footpoints were allowed (and actually had) to move poloidally. If, instead, one insisted 
on having the footpoints tied firmly to the disk surface, then the flux distribution on this surface 
(which, in the self-similar model, must be a power law extending to arbitrary large radii ) could 
not change, at least on the short, rotation-period time scale considered here. Then n and, hence, 
A9 would both stay finite. Under these circumstances the prospects for magnetic reconnection 
to occur in a timely manner would be very slim, at least in the two-dimensional framework (see 
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Uzdensky et al. 2002b). 

Thus, the transformation of the current concentration region into a true, infinitesimally-thin current 
sheet, even in the asymptotic sense, cannot, in general, be taken for granted. 

The failure of the self-similar force-free model to provide a plausible current-sheet formation (and 
hence reconnection) scenario makes it increasingly important to try to understand the asymptotic 
current-sheet formation process in a more general situation where the self-similar model does not 
apply. In particular, it is interesting to ask whether (still within a force-free framework) the current- 
concentration region can become very thin, thus indicating a way towards forming a true current 
sheet. If such a thin current layer does form, what is its structure? In particular, what is the 
^'-profile of the angular width and how does it depend on the twist angle profile A$(^')? 

It is important to realize that in the self-similar model there is no special radial scale and hence all 
the field lines must open up simultaneously. In a general, non-self-similar situation, however, there 
is a possibility of a partial field-line opening,^ with the domain of open field lines being adjacent 
to that of closed field lines. As we shall discuss below, this can greatly facilitate the current-sheet 
formation process. 

On the intuitive level, the basic idea of how a thin magnetic structure can form in the disk mag- 
netosphere is simple. Imagine a non-uniformly rotating disk. In particular, consider a system on 
the brink of the partial field- line opening: let there be a field line '^c such, that, as the system 
approaches a certain critical time tc, the field lines outside of *c (i-e., < *c) tend to open, 
while those inside stay closed."'' Now, if Ar2(^') [and hence the twist A<I>(^') = AQt] is nonuniform, 
say, rising outward, then the degree of expansion of field lines will also be non-uniform. Indeed, 
two neighboring field lines \l'i = ^'c -l- and ^'2 = ^ ^'^ will have their footpoints close to 
each other, while their apexes r'ap(^'i) and rap(^'2) will be very far apart since ^ap(^^'l) should stay 
finite whereas rap(*2) — 00 as t ^ tc- In other words, (ilograp(^')/cilog * at = ^'c will go to 
infinity as t ^ tc- As we shall show in this paper, a thinning of the apex region (i.e., current layer 
formation) is characteristic for a situation like this, which is consistent with the spirit of the studies 
by LBB94 and by Wolfson (1995) who both show that a current sheet forms asymptotically when 
n = d log ^'/d log r — > 0. 

Thus, the main thrust of this paper is to demonstrate how twisting of field lines leads to the 
formation of thin current structures that asymptotically become thinner and thinner as the field 
approaches a partially-open state in finite time. In § 2 we present the basic geometry of the problem 
and a description of our model. In §§ 3-5 we discuss the three main components of the model: 



*The importance of partial field-line opening in solar coronal processes has been recognized and emphasized by Low 
(1990) and by Wolfson & Low (1992), who suggested that a partially open field configuration may be energetically 
accessible oven when a completely open configuration is not. 

^We count poloidal magnetic flux ^' from the outside inward, i.e., we set ^(7r/2,r — > 00) — » 0. 
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§ 3 focuses on the purely geometrical relationship between the twist angle of the field lines and 
the relative strength of the toroidal field, while §§ 4 and 5 describe our treatment of the 0— and 
the radial force balance equations, respectively. In § 6 wc formulate and discuss the final set of 
differential equations describing the behavior of the angular thickness of the current concentration 
region. Also in that section we consider three specific examples: § 6.1 deals with the constant- twist 
case, § 6.2 investigates the behavior in the vicinity of the critical twist angle, and § 6.3 describes a 
numerical solution of our equations for the case of a locally-enhanced twist angle. We present our 
conclusions in § 7. 

Finally, we note that although we consider this problem in the accretion disk context, our methods 
and the main findings can be directly applied in the solar corona context. In particular, we note 
that the differential-rotation-driven process of partial field-line opening (and the detachment of 
the toroidal plasmoid associated with it) is essentially very similar to a coronal mass ejection; in 
addition, the partially opened field configuration considered in this paper is similar to a coronal 
helmet streamer configuration (see, for example. Low 2001). 



2. Description of the Model 

In spherical coordinates (r,9,(/)), the poloidal components of an axisymmetric magnetic field can 
be written in terms of the poloidal magnetic fiux function ^ as 

_ _J_dl ... 
' sine 89' ^ ' 

1 5* 

rsmt) or 

In a force-free equilibrium the magnetic flux function satisfies the Grad-Shafranov equation, 

sine d ( 1 , 



where the generating function 

F = B^r sine (4) 

is 2/c times the total poloidal current flowing through the circle deflned by G [0, 27r] at flxed 
{r,e). Since in an axisymmetric equilibrium poloidal current must follow poloidal field lines, F is 
constant along the field lines, i.e., F = F{'^). 

Now, consider a configuration with the dipole-like field topology, and assume that the magnetic 
axis coincides with the rotation axes of both the disk and the star. We are interested in a situation 
where the system is approaching the point of field-line opening, so that the field lines are expanded 
very strongly along the direction of their apexes e = e^p, so that rap(*) 3> ro(^) (where rap(^') and 
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Fig. 3. — Geometry of the problem: strongly expanded field lines in the magnetosphere of a 
magnetically-linked star-disk system. 

ro(^') are the radial positions of the apex and the disk footpoint of field line see Fig. 3). In this 
asymptotic regime 9ap(^) should not change much from one field line to another. In the self-similar 
model, for example, 0ap(*) is, of course, exactly the same for all field lines and varies by only a few 
degrees as n changes from to 1; and in the non-self-similar numerical model of Uzdensky et al. 
(2002a), 0ap(^) hes between 60 and 70°. In any case, the ^ dependence of ^ap is much weaker than 
that of the field line parameters rap and A9 (the latter is defined immediately below); therefore, 
we shall neglect this dependence in our analysis and assume that the apex direction is the same for 
all field lines. 



Next, for a given field line let us define the angular width A^i = A9{^) of the region enveloped 
by this field line. For definiteness, we define A9{^) in terms of the shape r{9, ^') of the field line as 



where e = 2.7183... is the base of the natural logarithm. 

In addition to A9i = A9{'^), we also define A92 as A6'2(^') = A0(*), where if is the field line with 

the apex radius r'ap(^') equal to e times rap(^) (sec Fig. 3). The quantity A^2(^) can also serve as 
an estimate for the characteristic angular width [at a given v — ^ap(^)] *^f the region where most 
of the toroidal flux, as well as most of the toroidal current, are concentrated. 

We are interested in the situation where thin structures are about to form in the magnetosphere. 
Therefore, throughout this paper we shall systematically assume that 



6'ap(*) = const. 



(5) 





A^(*) < 1, 



(7) 
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for all field lines * under consideration. 

Our main goal is to evaluate A^(^') and to determine how it changes in response to an increased 
twist angle. 

Our program has three main ingredients. First, there is a purely geometric constraint relating the 
twist angle A$(^') to the function and to the poloidal field Bq. The other two ingredients are 

dynamic — they refiect the force-free magnetic force balance in two directions, across and along 
the apex line (the 9 and the radial directions, respectively). Upon the completion of this program 
we shall obtain three equations relating three functions, F{^), A^(^'), and rap(^'), to each other. 
This will allow us to investigate under what circumstances thin structures can form. We start with 
the geometric constraint. 



3. The Twist Angle Constraint 

The twist angle A$('I') of a field line is defined as the difference between the toroidal positions of 
the footpoints of this field line on the surface of the disk and that of the star. It can be written as 

A$(*) = F{-^)I{-^). (8) 
Here, the function is given by an integral along the field line 



J Bersm^9 ^ ' 



Ber sin^ 9 



Using equation (2), we can express the poloidal field Bq in terms of the function r{^,9) as 

1 _ 1 f dr\ 

r sin 9 dr rsui9 \ J 



1 I fdr\~^ 

Be = -—r-,^ = (^.] . (10) 



Then, 

/(*) = - ~ - A 7^ 9) (11) 

^ ' J sm9 J sin^^ ' ^' ^ ' 

where we have neglected the dependence of 9* (the angular position of the footpoint on the surface 
of the star) on We can do this because, in the regime of interest to our present study, only a 
relatively small amount of twist resides near the disk and the star surfaces, while the dominant 
contribution to the integral /(^^') and, hence, to the total twist for a strongly expanded field 
line comes from the part of the field line near the apex, 9 = 9a,p. This is justified by the property 
of the twist to propagate along the flux tube and concentrate in the region of the weakest field (see 
Parker 1979), which in our case coincides with the apex region. 




apex 



apex 

e 







2Ae 



e 



ap 



e 



(b) 



Fig. 4. — The two examples of field-line shapes used in the text to demonstrate the calculation of 
the parameter k: (a) the triangular shape; (b) the rectangular shape. 



Furthermore, we can see that, if a field line is expanded strongly, then we can write 

it/2 



4^ri^,e)^KA9{^)'-^. (12) 
J sm0 sm6'ap 



This equation serves as the definition of the parameter k. Basically, this parameter describes the 
shape of the part of the field line near the apex [i.e., the part giving the main contribution to 
independent of the angular size A^(^) and the radial size rap(*), which appear in equation (12) 
explicitly. If all the field lines keep the same basic functional form near the apex, while only the 
typical angular and radial scales differ, we can regard k as being independent of ^. The validity 
of the assumption k = const is further strengthened by the observation that k varies very little 
over a broad range of field-line shapes. To demonstrate this, we here explicitly calculate, as an 
illustration, the values of k for two extremely different shapes, the triangular shape (Fig. 4a) and 
the rectangular shape (Fig. 4b). 

In the triangular case, shown in Figure 4a, we have, using the fact that A9 <C 1 and hence 
sin 9 ~ sin^ap = const. 



/ 



„/,T, m l'-ap(^) 



sin 9 2 sin ^ap 1 — 1 /e 

and so 



In the rectangular case, shown in Figure 4b, we again use <C 1 to approximate sin^ ~ sin^ap = 
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const, and obtain 



and so 



de 



suit 



rap(^) 

sin 6: 



AO, 



ap 



Kb = 1.0. 

Thus, from now on, we shall assume that k is constant. 
Using equations (11) and (12) we can now write 

K d 



sin ^ap 



(A^rap) . 



(13) 



Substituting this into equation (8), we get 

A$(*) sin^ap 



F(*) = - 



Krap(^') Ae{-^) {d\ogA9/d\og-^) - S(*) ' 



where 



d log rap (^') 



dlog* 



dr. 



ap 



(14) 
(15) 



4. The Force Balance 



Consider the 9 component of the force balance equation: 
47r 47r / \ 

— [j X B]e = — [j^Br - jrB^j 

= —-K-irBe) ^ — Ism^S,^) = 0. (16) 

r or r od rsmO od \ / 



Here, the first term represents the 6 component of the magnetic tension force, and the second and 
the third terms represent the 9 derivatives of the radial and the toroidal field pressures, respectively. 
Intuitively, it is clear that the first term must be small compared with the other two terms, because 
the primary force balance across the apex region should be between the toroidal field pressure at the 
apex and the radial field pressure outside the apex. The magnetic tension acts mainly in the radial 
direction while its 9 component is small. The same conclusion can also be reached upon noticing 
that the first two terms in equation (16) come from the corresponding two terms in the expression 
for the the toroidal current density, = {l/r){d/dr){rBg) — {l/r){dBr/d9). It is intuitively clear 
that, in the configuration where the field lines are strongly elongated in the radial direction and the 
radial field changes its sign across the apex line over a very narrow region, the second term gives a 
much greater contribution to than the first term. In § 6, we shall come back to this conclusion 
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and make a more rigorous a posteriori check of its validity, but now we shall just neglect the first 
term in equation (16) and thus get 



smOde 



r\ 7—) 

(^^sine) +5,-^ = 0. (17) 



Notice that i?</, changes on the angular scale A9 1, while sin^ changes on the angular scale of 
order 1 rad. Therefore, we can neglect the variation of sin 9 with 9 and write 



89 



+ =0 



Bi + Br = const 



at fixed r. 



(18) 



At the apex 9 = 9a,p, Br = 0, while outside the apex region, 9 > 9a,p + A^2/2, B^ becomes small. 
Hence, we obtain (consistent with the intuitive explanation given above) the pressure balance 
between the toroidal field at the apex and the radial field pressure outside the apex: 



Br 



ap) ' 



> 9^p + 



A9, 



(19) 



The toroidal field B^{ 

'"apj^ap) is easily expressed in terms of F(^): 



-B0(rap,^ap) 



F(*) 



rap(*)sin6'ap' 



(20) 



As for Brirg^p, 9 > ^ap + A02/2), we shall estimate it as follows. The total poloidal flux A^' through 
the interval 9 G [^ap + A^2/2,vr/2] per one radian in the toroidal direction at fixed r = rap(^^') 
is equal to A^' = *(r, ^ap + A^2/2) — *(r, 7r/2). But, since we assume that A^2 ^ 1, we can 
approximate *(rap(^'), ^ap + A02/2) ~ ^'(rap(^'), ^ap) = ^- In addition, for strongly expanded 
field lines, rap(*) ^ ro{^), and hence ^ S> ^'[ro = rap(*)] (where we take ^' = at infinity). 
Thus, we can evaluate A^* as A^* « ^. 

On the other hand. A* can be estimated as 

7r/2 

A* = -r2p(*) J Brsm9d9. (21) 

6»ap+A6»2/2 



In the interval under consideration (i.e., 9 G [^ap + ^^2/2, 7r/2]), Br is roughly uniform, and so, 
using A* we can write 

i?.(rap,^>^ap + ^)--^<0, (22) 
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where 

7r/2 

aj"^ = J sm9de = cos + ^) cos ^ap 2 (23) 

6»ap+Ae2/2 

is a constant. 

Note that estimate (22) is consistent with the self-similar model, which would give (see §§ 2.1 
and 2.3 of Uzdensky et al. 2002a) 

Bf"-'-'-(r.p,» > + ^) = ^g"'"/^' ^ . (24) 
^ 2 J Go,maxSm6'ap r^p(5') 

Since the function Go{9) is roughly linear in the interval [^ap + ^6*2/2, 7r/2], we can write Go,max — 
-G'o{Tr/2) (7r/2-^ap), and so 



Br (rap, ^ > ^ap + - ^ap ' T^i^) ' ^''^ 

Thus, combining equations (19), (20), and (22), we get 



Combining this with equation (14), we get 



F(*) = aisineap:r^-i-T- (26) 



dlogAg(^) A<&(^) 

dlog* ^ ^ KaiAe{-^)' ^ ' 

This equation presents a physically sound and adequate description of the system's behavior, with 
most of the uncertainties associated with its derivation being hidden in the constants k and oi . 



5. The Radial Force Balance 



Consider now the radial force balance along the apex line 9 = ^ap • Physically, we have here a balance 
between the radial outward pressure of the toroidal field, the term on the right hand side, and the 
poloidal field tension, second term on the left hand side of the Grad-Shafranov equation (3). The 
first term on the left hand side, the poloidal field pressure, is small in our limit <IC 1. The fact 
that this first term is much smaller than the second term can also be understood as the statement 
that, in the expression for the toroidal current density, = {\/r){d/dr){rBQ) — {l/r){dBr/d9), 
the first term is much smaller than the second term. Thus, we see that the criterion for validity of 
the assumption that the first term in equation (3) is much smaller than the second term coincides 
with the criterion for validity of a similar assumption we made when deriving equation (17) from 
equation (16). Thus, both these assumptions will be justified simultaneously in the next section. 
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So, we now have 



sin^ap a / 1 9* 



= -FF'{-^), . 



(28) 



apex 



(Note that our interpretation of the Grad-Shafranov equation only makes sense if the toroidal field 
pressure is directed outward, which corresponds to dF^/d"^ > 0). 



Neglecting, as usual, the 9 variation of sin 6 in comparison with that of the magnetic field quantities, 
we write 

1 d'^^ 

= -FF'{^),. (29) 

apex 



rip 89^ 



We can estimate d'^^/d9'^ at 9 = 9ap as 



d9^ 



a2 



apex 



A92/2 



r=ra.p 



where 02 is another constant of order one. 
At the apex d'^ /d9 = by definition, so 



apex 



2a2 

A9^~d9 



6'ap+A6l2/2 



■^r^p sin 6'apSr (rap,9 = 9ap + 



(30) 



(31) 



Using our estimate (22) for Brir^p, 9 = 9ap + A02/2), we get 



52^ 



89^ 



apex 



■^^sm^ap. 



(32) 



Then, from equation (29), 



FF'{^) 



2aia2 ^ 
A92 



ap 



(33) 



Thus we have obtained the following system of equations: 



F{^) 
FF'(^') 



dlog A9{^) 
dlog* 



ai sm C^ap 

2aia2 ^ 

A^>(^') 
~KaiA0(*) 



rap(*) ' 



(34) 
(35) 
(36) 



We view this system as a system of three equations for three functions, F{^), rap(^'), and A9{^). 



-15- 



To simplify this system, we shall first use equation (34) to eliminate function F{'^). We have 



(37) 



Substituting this, together with equation (34), into equation (35), we get 

2a2 



l + 5(*) A^2(*) 



= const = 0(1) . 



ai sm- (^ap 

For example, for ai = 1/2, a2 = 1, and sin^ap — 1, we have C ~ 1/4. 

Now note that, because we are in the small-A^ regime, from equation (38) it follows that 

f^log rap 



(38) 



5(*) = 



> 1. 



Then equation (38) becomes 



dlog* 
C5(*)A^2(*) = 1 • 



(39) 
(40) 



It is interesting to note that, since S = — dlograp/dlog then ^(r, ^?ap) ^ i" ^^^i that is, 1/5" 
represents a direct analog of the power exponent n in the self-similar models of VB94, Uzdensky 
et al. 2002a, Low & Lou (1990), and Wolfson (1995) (and the power exponent p in LBB94). Thus, 
the last equation can be interpreted as a statement that n= 1/S ^ A62 — > as a thin current layer 
forms, in agreement with LBB94 and Wolfson (1995). This shows that a current sheet formation 
is intrinsically related to the increase in S, and hence can be attributed to a growing difference in 
the expansion ratios of neighboring field lines with different twist angles, as discussed at the end 
of 5 1. 



6. Analysis and Interpretation of the Results 

Thus we managed to reduce the problem to a system of two equations, namely, equation (40) and 
equation (36). The latter can be rewritten as 

1 dA0 A0 . , 

= -u{x). (41) 



SAe2 dx A02 
Here 

x = log^', (42) 

and 

u{x) ^ (43) 



-16- 



where 



~c' 



sm 



const = 0(1) 



(44) 



is some critical twist angle. For example, for k = 1, a2 = 1, and sin^ap — 1) we get A$c ~ 2 rad. 
As we shall see below, the behavior of the system below and above A$c greatly differs. 



It is convenient to define a new independent variable 

t = lograp(^'). 

Then, 



and equation (41) becomes 



Aeit) 

A9{t+1) 



S = —dt/dx, 
^ dlogM 



dt 



u[xit)] 



(45) 
(46) 

(47) 



At this point we would like to pause in order to verify the validity of two steps in our program, 
namely the derivation of equations (17) and (28) in § 4 and § 5, respectively. The validity of 
both these steps was based on a single assumption that, at the apex 6 = ^ap the first term in the 
expression for the toroidal current density, j^f, = {l/r){d/dr){rBg) — {l/r){dBr/d6), is much smaller 
than the second term (in the limit A6 <C 1 that we consider here). At this stage in our analysis 
we now have at our disposal all the tools necessary to estimate the magnitude of these two terms. 
Indeed, using equations (2), (39), and (40), the size of the first term can be evaluated as 











apex 


Qj.2 



apex 



^Ip(*) 



o 



+ 



dA9<2 



dt 



' ap 



(48) 



At the same time, using equations (1), and (32), the second term, {l/r){dBj./dO), can be estimated 

(49) 



at ^ = ^ap as 



1 dBr 


1 




r 86 


— 2 

apex ^ 


09^ 



1 



T 

apex ap 



> 



Thus, we see that the approximation we used is indeed well justified in the small-A6' limit. 
We now turn to the application of our results to several particular situations. 



6.1. The Constant-u Solution 

The non-local nature of equation (47) makes it very difficult to analyze. The situation is complicated 
even further by the fact that the dependence x{t) and, hence, the right hand side of equation (47) 
are not explicitly known. However, an exact solution of this equation can be found in the case of 
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a disk rotating as a solid body, A$(*) = const. This may be, in fact, a good approximation for 
the outer parts of a Keplerian disk, r ^ (where is the corotation radius), where Ar2(^') = 
^disk(^) — ~ — f^* = const. In this uniformly-rotating disk model, one thus has u(x) = const. 
The solution that satisfies the boundary condition A.6{to) = AOq is 

A0(t) = A0oe(^"«^(*"*°\ (50) 

and the constant ^ is uniquely determined by u via 

n = ^e(«-i\ (51) 



Next, using equation (40), we find S = S'o exp[(^ - l)(i - io)], where = S{to) = u/C^AOo. Then, 
using equation (46), we can calculate x{t), 

x = xo + ^-^(e(i-«(*-*'')-l), (52) 

and thus determine AO as a function of the poloidal flux: 

u ^ — 1 ^ 1 ^' 

A9ix) = AOo + p^(^ - l)(x - xo) = AOo + ^e^'' log ^. (53) 



For each u there is a unique solution S^, such that ^{u > 1) > 1 and ^{u < 1) < 1. This means that, 
if M > 1 (that is, A$ > A$c), then AO decreases outward, and if w < 1 (A$ > A^>c), then AO 
increases outward. The special case u = 1 corresponds to ^ = 1 and A^ = const and thus describes 
the approach to the singularity A$c in the self-similar case discussed by VB94 and Uzdensky et 
al. (2002a). 

Notice that in the case n < 1 (^ < 1), AO increases outward (with decreasing *) until our model's 
main assumption AO <^ 1 breaks down. This condition defines the outer boundary ^'min of appli- 
cability of the present model, which can be estimated roughly as the value of for which AO ~ 1: 

(54) 

We see that the range of validity of our model widens (^'min decreases) as ^ — > 1. 
In the opposite case n > 1 (^ > 1) there exists a certain critical field line, 

(55) 

on which AO vanishes while rap becomes equal to infinity. If u (and hence ^) just barely exceeds 1, 
then this point is very far away: *c ^ and r-o(*c) ^ ?'o(*o)- As u and ^ increase, also 
increases and the critical field line moves inward. 



*min(0 - *oexp 



C 



= ^'o exp 



CAO, 
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We interpret the vanishing of A6{'^c) as a partial field-hne opening, with ^'c marking the boundary 
between open and closed field lines. Our model describes only closed field lines and thus is valid 
only for ^ > For any given field line ^ > '''ap(^) is still finite, and, since rap(^'c) = oo, the 
field line ^ with rap(^) = erap(^') still lies between ^ and ^'c- We can see then that the analysis 
given in the previous sections will still be valid for any field line * > no matter how close 
to '^c- For open field lines, * < ^c, our treatment is, of course, no longer valid. However, such an 
elaborate analysis is not needed for these field lines because their physics is very simple: the open 
field < *c) is potential, = 0, and mostly radial. 

We would like to remark that expression (53) describes the behavior of the angular width of the 
current concentration region as a function of We can also see whether the physical width of 
this region (which probably is of greater interest, as far as reconnection is concerned) decreases 
as —>■ 0. Prom equation (52) we have 

= [Soi^ - l)(x - xo) + 1]^ , (56) 

and so 

rap(*)Ae(*) DC [Ae{^)]^i. (57) 

Wc sec that TapAfl with A9 if ^ > 2, i.e., if n > 2e ~ 5.4. Thus, wc envision the following 
scenario. As the twist is increased beyond u = 1, there will be a line beyond which the field will 
be open: T'ap(^ — ^ ^c) — ^ -^s the twist is increased even further, more and more flux becomes 
open. At the same time, as u is increased beyond 2e, then, according to the power law (57), not 
only the angular, but also the actual physical width of the current concentration region, rapA0, 
shrinks to zero as one takes the limit ^' This suggests that increasing the twist angle well 

above the critical value can not only result in a greater amount of open fiux, but also can lead to 
a narrower current layer, thus facilitating magnetic reconnection. 

We would like to emphasize that, while our analysis predicts a finite-time partial field opening, it 
does not, strictly speaking, predict a finite-time current-sheet formation in a true sense. Indeed, at 
any given > A$c, and for any given still-closed field line ^ > ^c(A^*), the current concentration 
region has a non-zero angular thickness A6{^,t) (which, however, tends to zero as * — > ^c)- If 
one fixes a field line and monitors how the angular thickness A^(*,t) changes with time, then 
one finds that it does indeed go to zero in finite time at the moment of opening of this field line. If, 
however, one looks at any given finite radius r (instead of fixing a field line then the situation 
is very different. Since rap(^'c) = oo, there is always a closed field line ^' > with rap(^') = r; 
thus, the angular thickness A9{r, t) stays finite at the moment of first field-line opening (when 
A$ = A$c)- As the twisting continues, A9{r, t) will probably decrease and may in fact go to zero 
asymptotically as t ^ oo. Thus, one may say that even though one has a finite-time partial field- 
line opening, there is no true finite-time current sheet formation. This is indeed a very important 
distinction. Of course, on a practical note, what is important here is the fact that the current layer 
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becomes very thin: at some point during the thinning process finite resistivity may come into play, 
leading to field-line reconnection. 

Finally, we need to point out that, when partial field opening takes place, some parts of an opening 
field line move very rapidly (e.g., the velocity of the apex approaches infinity at the point of 
opening). We have to make it clear that these infinite velocities are an artifact of the equilibrium 
assumption. They should be understood only in the sense that they become very large compared 
with the footpoint rotation velocity. In reality, the expansion velocity will be limited, most likely, 
by finite inertia of the plasma, i.e., by the Alfven speed. Another concern is the assumption that the 
suggested sequence of partially-open field configurations is continuous. It is actually not completely 
clear whether this always has to be so. In reality, it may be possible that, upon reaching a critical 
twist, a quasi-static evolution along a sequence of force-free fields will end abruptly by opening up 
a finite portion of the fiux. 



6.2. Solution in the vicinity of A$c 

Another important advance in our understanding of equation (47) can be made in a situation where 
A<I>(^') is nonuniform and passes through the critical twist angle A$c (for example, as one moves 
outward, i.e., as ^ is decreased). 

Let us consider a close vicinity of the point xq = log^o where u{xo) = 1. In this region we can 
Taylor-expand u{x) and keep only the linear term: 

du 

u{x) = 1 — r]{x — xo), where r] = ^^(^o) > 0- (58) 
Let us define two functions, 

and 

(60) 

and consider a region around to = t{xo) corresponding to a close vicinity of xq, such that |«— 1| <^ 1 
within this region. We anticipate that varies relatively little so that \e{t)\ ^ 1 everywhere in 
this region. 

Then equation (47) becomes 

{l-e){l-a) = u^l, (61) 
and therefore \cr{t)\ <IC 1, which gives us a relationship between e{t) and cr{t): 



e{t)+a{t) = l-u{t) < 1. 



(62) 
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A second relationship follows from the definitions of e{t) and cr(t): 



t+i 

z{t) ~ J a{i)dt, 



(63) 



where we have neglected higher-order corrections in a{t). Using this expression we can differentiate 
equation (62) and get 

a'{t)+a{t + l) -a{t) = -u'{t), (64) 
where prime denotes the derivative with respect to t. 

Now let us Taylor-expand a{t + 1) around the point t and keep only the first two terms: 

a{t -Fl) « a{t) + a'{t) . (65) 

Then we can write 



a'{t) 



u'{t) 



(66) 



Wc now need to make a little digression in order to discuss the validity of this step, that is of 
keeping only the first two terms in the expansion of a{t + 1). In particular, we need to verify 
that we can safely neglect the next term in the Taylor expansion, 0.5 a" (t). This term needs to be 
compared with the last term that we kept in the expansion, that is with (7'{t), which, according to 
equation (66), can be rewritten as 



u'it) 



1 du{x) 
2S 



C du{x) 



Ae{t + 1) , 



2 2S dx 2 dx 
where we made use of equations (46) and (40). At the same time, the first neglected term is 



(67) 



a"{t) d'^u{x 



C du(x) 



dx^ 



4 dx 



where we again used equations (40) and (46), as well as equation (60). The second term in this 
expression is manifestly much smaller than (67) because we here are considering a region where 
\a{t)\ <IC 1. Thus we see that our approximation (65) is valid if 



du{x) 



dx 



rj > 



d^u{x) 



dx"^ 



M{t + l). 



(68) 



The function u{x) is the prescribed twist angle as a function of flux; generally, it varies on a scale 
that is finite in ^ and, hence, in x. This means that both dujdx and (Pu/dx"^ are of order 1; 
therefore, since we work in the regime AO <C 1, the above condition is easily satisfied and thus our 
truncation of the Taylor expansion (65) is valid. 



We now continue with our analysis. Since we consider a region where is almost constant, we can, 
to the lowest order in cr, substitute A0{t + 1) in equation (67) by A^o = ^^(^o)) where to is defined 
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by x{to) = Xq. In a similar vain, du/dx also varies very little in the region under consideration, so, 
using equation (58) we can substitute it by —77. Then we can integrate equation (67) to get 



a{t)=ao~^Mo{t-to) + 0{a^), 



(69) 



where ctq = cr(to)- 

Combining this expression with the definition of cr(t), we get the following expression for AO(t): 



Ae{t) 



A9o exp 



croit - to) 



■ exp 



-AOoit - to)' 



~ Aeo[l + ao{t-to)- ^AOoit-to)'^ 



(70) 



Then, from equation (63) it follows that 



Cri 



and from equation (62) we get 



e{t) = ao-^Aeo[l + 2{t-to)], 



Cn 

(Jo = a{to) = -^A9o > 



(71) 
(72) 



Thus, we can write 



Ae{t) = Aeo 1 + ao{t - to)[l - 2{t - to)] 



(73) 



We see that Ad{t) behaves regularly at t = to; for t > to, it first increases, reaches a maximum at 
t = to + 1/4 (i-c, X — Xq = —CA60/4 <C 1), and then starts to decrease. If u > 1 for all x < xq, 
then, guided by our constant-u solution, we can predict that at some point will reach zero, 
implying a partial field-line opening. 



6.3. Locally enhanced twist 



As a last example, we present the results of numerical solution of the system (47)-(46) for the case 
of a locally enhanced u{x). 

We consider an interval x G [xout ) a^in] corresponding to G [*out ) ^in] • In order to be able to make 

a connection with the constant-u analytical solution, we take the function u{x) to consist of three 

pieces (see Figure 5): 

u{x) = uo <l, Xout <x<xi; 

u{x) = Ul > l,Xl < X < X2; 

and 

u{x) = Uo,X2 < X < Xin. 
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u(x) 



out 



Fig. 5. — The u{x) profile investigated in § 6.3. 



Here, uq = const < 1 is some background value, which we keep fixed, and Au = ui — uq = const 
is the localized enhancement. We study a family of solutions of the system (47)-(46) with a fixed 
value of uq, uq = 0.655.. [which corresponds to = C(^o) = 0.8], but with several different values 
of 71,1. We start with ui = 1.2, and then gradually increase it until the sequence terminates at 
ui ~ 1.8. In order for us to be able to make a meaningful comparison of the solutions with different 
values of ui, we fix the value of AO on the innermost field line, A9{xm) = AOq. 

Now, based on our experience with constant-n solutions, we can naively expect the resulting function 
A6{x) to be roughly a piece- wise linear function, with A6{x) decreasing linearly in the intervals 
Xout < X < xi and X2 < x < Xm [with the slope uo(^o — l)/C'^o], and increasing linearly in the 
interval xi < x < X2 [with the slope ^1(^1 — 1)/C^i, where ^1 = C{ui) > 1]. We thus expect AO to 
have a maximum at X2, A^^^ = A0{x2) = AOq + uo{$,o — 1.){x2 — Xi^/C^Q, and a minimum at xi, 
Ae(xi) = AOt^ - - 1){X2 - xi)/Ca. 



mm 



Note that, when u\ reaches a certain value ui^max 
such that 

(6,max - 1) exp(^i 

,1 



> 1 [and, correspondingly, ^i,max = C(^ti,max) > 1] 

„a.-l) = ^^, (74) 



X2 - Xi 

then A^^*^ becomes zero. We interpret this as the point in the sequence of equilibria at which a 
partial field-line opening occurs at = Note that ui^max, i-e., the maximum allowable value 
of Ml, depends inversely on the size {x\ — X2) of the enhanced- twist region. We thus see that, if 
only a portion of field lines (those between and ^2) are participating in twisting and expansion, 
then the thinning of the current layer, and the corresponding partial field-line opening, is achieved 
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not at -u = 1, but at a somewhat higher value. 

When performing the actual computation, we advance along t from tout (corresponding to the 
outermost field line *out) inward. We first pick a guess for the initial value of at the outer 
boundary, Agouti and then use our const-'u solution to prescribe the function A9{t) and hence x{t) 
in the initial interval of unit length, t G [tout ~ Ij^out]- After that, we proceed to integrate the two 
differential equations: 

^=A»(t)-u[i(t)]A«(t+l). (75) 



We stop when the innermost field line of the domain, x = Xi^, ^ = ^m, is reached and compare the 
computed value of A^ at this point with the prescribed value A9i^ = A9q. Wc then iterate with 
respect to Agouti i-e., we change A^out and repeat the procedure until A0(xin) finds itself within a 
small vicinity (which we usually take to be 1%) of the desired value A^o- 

In Figure 6 we present the results of our computations for A^o = 0.3, Xout = 0-0) xi = 0.5, X2 = 1.5, 
Xin = 2.0, C = I, uq = 0.655, and ui between 1.2 and 1.8. Wc indeed see that the obtained solutions 
A9(x) are in good agreement with the analytically-predicted picture described above. In particular, 
we find A^max — 0.43, and iii,max — 1-8- This is very close to the naive prediction (74), which 
gives A^max = 0.382, ^i,max = 1-287, and uf^^ = 1.71 for the set of parameters that we chose. 
And, if we use the actual value A^max = 0.43 in equation (74), we then get ^i,max = 1-314, and 

^!%ax = 1-80- 

At the same time we observe that, as ui^^^x is approached, the apexes of the field lines move out to 
infinity explosively fast, signaling a partial field-line opening at a finite twist- This opening can be 
illustrated by the behavior of the natural logarithm of the ratio of rap of the outermost field line to 
Tap of the innermost field line, i.e., by tout — t^i = log[rap(^'out)/''ap(^m)]- Figure 7 shows this ratio 
as a function of ui. We thus see an unambiguous evidence that the thinning of the current layer 
occurs simultaneously with the partial field-line opening, in accord with our expectations discussed 
in 5 1. 



7. Conclusions 



In this paper we addressed the issue of the formation of thin current structures in the magneto- 
sphere of an accretion disk in response to the field-line twisting. We presented a simple analytical 
model that illustrates the asymptotic behavior of a force-free axisymmetric magnetic field above 
a thin conducting disk in the limit when such structures are about to form. Wc showed that 
there exist a finite critical twist angle, A^c, such that the expansion of the poloidal field lines 
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Fig. 6. — Angular thickness A9 as a function of x = log* for a range of values of ui between 
ui = 1.2 and ui = 1.8. 





Fig. 7. — Natural logarithm of the ratio of rap of the outermost field line to rap of the innermost 
field line, iout ~ ^in = log[rap(*out)/?^ap(*in)], as a function of the twist ui. 
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accelerates dramatically as A$c is exceeded. At the same time, in the case of a locally-enhanced 
twist, ^ a radially-extended thin conical current layer forms. Simultaneously, the field configuration 
approaches a partially-open state in finite time, with the apexes of a portion of the field lines going 
out to infinity. 

The presence of a thin current sheet is usually considered a pre-requisite for reconnection of magnetic 
field lines. It is also believed that reconnection will occur whenever a thin current layer is formed. 
We therefore believe that magnetic field lines that have become effectively open as a result of their 
twist-driven inflation may subsequently reconnect via the thin current layer that has been formed 
simultaneously with their opening-up. This conclusion leads to important consequences regarding 
the time evolution of magnetospheres of magnetically-linked star-disk systems. In particular, we 
believe that, once the magnetic link between the star and the disk is broken by the partial field- 
line opening, a current layer will form along the separatrix between the stellar field lines and the 
disk field lines and, as a result, the link will be re-established through reconnection of these field 
lines. Subsequently, if reconnection process is sufficiently fast, the field lines will contract to a 
less-stressed configuration allowing a new cycle of twisting, inflation, opening, and reconnection 
to begin, as was suggested by Aly & Kuijpcrs (1990) and by VB94 and illustrated in numerical 
simulations by Goodson et al. (1999). This quasi-periodic scenario is characterized by very rich 
physics with non-steady and violent behavior, perhaps not too different from that of the solar 
corona (e.g.. Low 2001). It can provide avenues for understanding such phenomena as disk winds, 
time- variable, knotted jets, episodes of rapid accretion, and variable accretion torque on the central 
star. 

One should note, however, that a true opening will be preceded by the plasma inertia becoming 
important. The inertial effects will tend to retard the expansion, since the pressure of the toroidal 
field pressure will have to work not only against the poloidal field tension, but also against the 
plasma inertia. Hence the inertial effects will in effect act against the toroidal field removal and 
will not help to form a current sheet. At some point during the expansion they will need to be 
taken into account and, even in the absence of reconnection, a true finite-time field opening will 
not happen. Instead, one will have a transition from the force-free regime into the inertial regime 
(wind regime), which would require solving the full set of MHD equations, although probably under 
some simplifying assumptions. 

I am very grateful to Stanislav Boldyrev, Arieh Konigl, and Bob Rosner for very useful discussions 
and comments. I am also grateful to the anonymous referee for his or her deep and insightful 
criticism that was very beneficial for the paper. I would also like to acknowledge the support by 
the AS CI/ Alliances Center for Astrophysical Thermonuclear Flashes at the University of Chicago 



An important example of the locally-onhanccd-twist case is a Keplerian disk whose inner edge is significantly 
closer to the central star than the corotation radius. In particular, if rinner < 0.63 Tco, then the relative star-disk 
angular velocity, Af2(r) = f2diBk(r) — f2*, reaches higher values in the localized inner region (rinner < r < rco) than at 
large distances (r > rco)- 
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